Coulomb scattering cross-section in a 2D electron gas and production of entangled 

electrons 
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Abstract 

We calculate the Coulomb scattering amplitude for two electrons injected with opposite momenta in an interacting 2DEG. 
We include the effect of the Fermi liquid background by solving the 2D Bethe-Salpeter equation for the two-particle Green 
function vertex, in the ladder and random phase approximations. This result is used to discuss the feasibility of producing 
spin EPR pairs in a 2DEG by collecting electrons emerging from collisions at a n/2 scattering angle, where only the entangled 
spin-singlets avoid the destructive interference resulting from quantum indistinguishability. Furthermore, we study the effective 
2D electron-electron interaction due to the exchange of virtual acoustic and optical phonons, and compare it to the Coulomb 
interaction. Finally, we show that the 2D Kohn-Luttinger pairing instability for the scattering electrons is negligible in a GaAs 
2DEG. 



I. INTRODUCTION 

Recent experimentsi*^ have allowed for the imaging of 
the coherent electron flow in a two-dimensional electron 
gas (2DEG), demonstrating a roughly directional injec- 
tion through a quantum point contact (QPC) tuned to 
his lowest transversal mode. We propose here to use 
such a setup to investigate Coulomb scattering in 2D, 
by measuring the scattering cross-section. This provides 
a natural motivation for solving a long-lasting problem 
in Fermi-liquid theory: finding the scattering amplitude 
for the Coulomb interaction in a 2D system, including 
the effect of the many-particle background of the inter- 
acting Fermi sea. We derive^ the scattering amplitude 
/ by solving the Bethe-Salpeter equation in the ladder 
approximation and for electrons in the Cooper channel 
(opposite momenta)^. This solution provides a use- 
ful addition to Fermi-liquid theory appHed to electron- 
electron interaction. The development of Fermi-liquid 
theory, which goes back over many decades, includes dis- 
cussions of screening^, the lifetime of quasi-particlesiA, 
the renormalization factor Z of the Green function2*ift, 
the effective massiSiii, and scatteringiSiiii4. An addi- 
tional issue is how strongly the Coulomb scattering is 
affected by the effective electron-electron interaction me- 
diated by the exchange of virtual phonons, which have 
been studied, for instance, in the context of Coulomb 
drag (see e.g. Refii^) and screeningi^. Another extension 
concerns the effect of lower dimensionalitjii^iii in Kohn- 
Luttinger superconductivitj^, and the question of the 
strength of superconducting fluctuations (if any), which 
could, in principle, spoil the Coulomb scattering. 

A second motivation for this work comes from the 
current efforts devoted to soHd-state implementations of 
quantum information protocols using the spin of indi- 
vidual electrons as qubitsiS'ifi'^. In particular, the ex- 
perimental demonstration of entangled (EPR) pairs of 
spin-qubit is still a present-day challenge, and has mo- 



tivated a number of theoretical proposals for entanglers, 
i.e. devices creating mobile (spin-) entangled electrons. 
These proposals relied on energy flltering via quantum 
dotsiSiSiiSSiSi and carbon nanotubes^, and/or the use of 
superconductor82iiSiiSSi2£. Other schemes are in closer 
relation with optics, and use beam-splitters for spinal, 
orbital^ or particle-hole entanglement^Si^. 

We proposal here a simple idea for the creation of 
spin-entangled pairs in a two-dimensional electron gas 
(2DEG), inspired by the well-known interference effect 
found in vacuum for the scattering of indistinguishable 
particleSii^. Using the fact that electron pairs in the 
singlet (triplet) spin state behave like spinless bosons 
(fermions) in spin- independent collisions, we propose to 
collect electrons emerging from electron-electron collision 
with a scattering angle 9 — 7r/2. In this situation, the 
destructive interference is complete for triplets; hence the 
collected electrons must be in the entangled singlet state 
15-) = (I fl) - I |t))/\/2, which is one the EPR states 
desirable for quantum information protocols. The ques- 
tion now arises whether this two-particle exchange effect 
survives in the presence of a sea of interacting electrons. 
Using our solution for the scattering amplitude /, we 
will show that the entanglement created (or rather, post- 
selected) by the collision should be observable in a real- 
istic 2DEG. 

We emphasize that our study of electron-electron in- 
teraction is motivated by the prospect of seeing experi- 
mental control on the propagation of electrons in 2DEGs 
and on their quasi-particle properties. The experiments 
of Refpii^ have used electron scattering off an SPM tip to 
image both quantum interferences and ray-like propaga- 
tion of electrons, including diffraction. They have also 
shown control over the quasi-particle lifetimes for hot 
electrons, in good agreement with Fermi-liquid theory. 
The theoretical work presented here indicate that such 
experimental studies can be extended to more general 
Fermi-liquid effects involving two quasi-particles^^. 
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Figure 1: Setup. Two quantum point contacts (QPC) allows the injection of electrons from two reservoirs with initial momenta pi ~ — P2- 
(a) Measuring the conductance as a function of the SPM tip position gives an estimate of the electron fiu5^ and, thereby, of the differential 
scattering length X{9) (top). Alternatively, one can define "bins" spanning different angles, and collect the current in the drain contacts 
(bottom), (b) EPR setup. The electrons are collected in two detectors (with an aperture angle 2S6) placed such that only electrons 
emerging from collisions with a scattering angle around tt/2 are detected. Because of antisymmetrization, the scattering amplitude 
identically vanishes for the spin-triplet states, allowing only the spin-entangled singlets (EPR pairs) to be collected, (c) Scattering 
parameters. The initial (pi,P2) and final (pi,P2) momenta are connected by a circle of radius p' = p due to energy and momentum 
conservation, where the relative momenta are p = (pi — p2)/2,p' = (pi — P2)/2 and 9 = Z(p,p') is the scattering angle. The 
Cooper channel is defined by 2a — Z(pi, — P2) 0. 



We start in Section ^ by describing the envisioned 
setup and the mechanism for the production of EPR 
pairs. We write down in Sec. IIIII the problem in a 
Fermi Uquid approach, solve the Bethe-Salpeter Eq. in 
the ladder approximation, using RPA and considering 
the Cooper channel (opposite incident momenta). The 
solution for the scattering amplitude (the t— matrix) is 
written in Eq. l(38jl in terms of a Fourier series, with ex- 
plicit expressions for the coefficients. We study in Sec. 
IIVI the scattering cross-section, and address in more de- 
tails the issue of the production and detection of the EPR 
pairs. We investigate in Sec. 0the electron-electron in- 
teraction mediated by phonons, and show that it does 
not have a significant effect on the scattering. In Sec. 
IVII we show that no superconducting instabilities arises 
from the Kohn-Luttinger mechanism.^-. We finally con- 
sider in more detail the case of small in the Appendix. 
We derive in Appendix fXl analytical expressions for the 
scattering amplitude and its derivative at = ■7r/2. In 
Appendix El we develop a different calculation vaHd for 
very small r^, which is needed to estimate the contribu- 
tion of forward scattering states. 

II. SETUP AND PRODUCTION OF EPR PAIRS 

The setup for the study of Coulomb collisions in a 
2DEG is described in Fig. ^ Two quantum point con- 
tacts (QPC) tuned in their lowest transversal mode filter 
electrons escaping from two thermal reservoirs^ and al- 
low them to collide with incident momenta pi ~ — P2 
and final momenta , P2 . A way to measure the scat- 
tering length is to use an SPM tip and to record the con- 
ductance across the sample, which provides an estimate 
for the local electron fiuJ«; alternatively, one can define 



'TdIus" spanning different angles, and collect the current 
in the drain contacts; see Fig. ^a). For the production of 
EPR pairs, described in Fig. db), the electrons are col- 
lected at two detectors placed so that only coUisions with 
a scattering angle 9 within a small window 59 around 7r/2 
are collected: 9 £ [7r/2 — 59, t: 12 + 59]. We introduce the 
conserved total momentum P = pi + p2 = pi + P2, the 
relative momenta p = ^(pi — P2), p' = ^(pi — P2), and 
the scattering angle 9 — Z(p',p); see Fig. QJc). The 
most favorable arrangement is the Cooper channel^ 

P2 ^ -Pi, (1) 

as is yields conservation of the individual energies: pi ~ 
P2 — p'l — P2j where pi = |pi|. As a consequence, the 
scattering angle 9 ~ Z(pi, pi) can be easily determined, 
while the EPR pairs have the same energy and should 
therefore arrive in the detectors at the same time. We 
consider incoming electrons with small excitation ener- 
gies = ft^pf/2m — Ep <C Ep above the Fermi energy 
Ep = h'^kj,/2m of the 2DEG (m is the effective mass>24. 

Now we describe in more details the production of spin- 
entangled electrons. First, we use the fact that the two- 
particle interference is totally destructive for fermions 
colliding with a scattering angle 9 = n/2, while, in con- 
trast, the scattering of bosons is enhanced compared to 
the classical value. This is seen in the corresponding 
cross sections ctb/f — l/l^*) i /("■ ~ Second, the 

fermionic character of a pair of particles also depends on 
its spin-stateSiSii^: a spin-singlet electron pair 

1^) = (I Ti> - I iT»/V2 (2) 

behaves, in a spin-independent coUision, Hke a bosonic 
pair because of its symmetrical orbital wavefunction, 
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while the triplets 

|To) = (I Ti) + I iT))/\/2 ; |T±> = I TT>, I ii> (3) 

behave like fermions. It is then clear that a 7r/2 scattering 
experiment could distill the singlet part of uncorrelated 
pairs of electrons -at least in free space. 

It might seem surprising to be able to collect spin sin- 
glet states when starting from two electrons having no 
spin correlations (they come from two independent un- 
polarized reservoirs). However, a spin-singlet component 
is always present, as seen from the change of basis for the 
density matrix describing the two spin state: 



relation^^ 
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\aG ((T(T 



(4) 



where \aa') corresponds to the two-electron state where 
the electron injected from the first (second) reservoir has 
spin cr (cr'). 

A real detector has a small, but finite aperture an- 
gle 250 around Q — tt/2, so that triplets will always be 
present. To examine the efficiency of this collision entan- 
gler, we wih define in Sec. I^the ratio TZ between the 
number of scattered triplets and singlets, Nt/s- We will 
find 



Ns 



/'(V2) 



(5) 



and show that 1/7/1 1 at 6* = n/2. Therefore, the 
number of triplets (which we want to avoid as they can 
be in an unentangled product state |T±)) is neghgible for 
small 69. 



III. CALCULATION OF THE SCATTERING AM- 
PLITUDE 

We study in this section the scattering between two 
electrons that are both above the Fermi surface — as op- 
posed to the standard calculation of the electron lifetime 
due to scattering of one electron above the surface with 
all the electrons present below the surface^. We con- 
sider a clean 2D Fermi hquid with Coulomb interaction, 
neglecting impurity scattering (the mean free path can 
be around 10 fim, which is larger than the size L ~ 1 /im 
of the envisioned setupi) . The effect of phonons wih be 
considered in Sec. 



A. RPA and Bethe-Salpeter equation 

In 2D, the scattering amplitude / for two particles with 
I, relative momentum p is linked to the t— matrix via the 



m 



(6) 



with the corresponding scattering cross-section ("length" 
in 2D) X{e) = 1/(6*)^. The 2D Coulomb interaction in 



Vc(q) = j dre' 



2-Kel 



(7) 



where Cg = e^/47reoer. We included the dielectric con- 
stant Cr of GaAs and set p = kp for future comparisons 
with the scattering in GaAs. The corresponding exact 
scattering i-matrix is given by^ 



tc{e) 



sin|6'/2| TO r(l-i<r) 



J-iT/i~2i<; In | sine/2| 



(8) 

with r(a;) the Gamma function and — mef^/kpfi^ ■ This 
yields 



\tcm 



\J Txc, tanh(7r<^) 
msin|6'/2| 



(9) 



Note that \tc\ ^ Vc, contrary to the situation in 3D. 

Next we include the effect of the many-particle back- 
ground, and calculate the i-matrix in the presence of the 
Fermi sea. The i-matrix is given by the vertex func- 
tion t — r(ijJi,2 £.1,2) appearing in the two-particle 
Green function'*; see Fig. Ela). We note that the arrange- 
ment P2 = — Pi corresponds to the well-known Cooper 
channel, discussed for instance in the context of Cooper 
instabilitji^. Consequently, we adopt the approach used 
by Kohn and Luttinger in their work on intrinsic super- 
conductivity in a 3D Fermi liquidAi. However, our calcu- 
lation differs in two ways. First, we consider a 2D system 
where the screened potential is non-analytic (because of 
the modulus <? = |q| instead of in 3D). Secondly, we are 
interested in the scattering amplitude, while Kohn and 
Luttinger focused on the instability in the vertex aris- 
ing from spherical harmonics of the crossed ('exchange') 
diagram (see A3 below). 

The two-particle Green function in real space is 

G,..,,,.,(l',2';l,2) = (-*)'(rci,^iC2v^4,4. 

(10) 

with the notation i = (xi,ii), ai =t, i are the spin in- 
dices, and T the time-ordering operator. For a spin- 
independent Hamiltonian we can writ 0^1*^ 

G,i4,,,,(l',2';l,2) = G(1)G(2)(2^)3 

X [,5(1' - l)d,'^aM-, - ^(1' - 2)5,^,, 5,^,,,] 
+ (*A)G(l')G(2')r_(l',2';l,2)G(l)G(2), (11) 

with the single-particle Green function G(i) and 
the (anti-) symmetrized vertex r±(l', 2'; 1, 2) = 
r(l',2';l,2) ± r(l',2';2,l) expressed in terms of 
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the unsymmetrized T. We consider the singlet /triplet 
basis {IS*), |To), |T±)} and introduce the corresponding 
creation operators 



^S/To 



(1,2) - l/%/2(, 



C-| -J- Crj I -p Cl I C. 



(1,2) 



u^2T y ' 



(12) 
(13) 



which we use to define the related two-particle Green 
functions 



G5/To(l',2';l,2) 



'ra5/To(l',2')4/^„(l,2)^ 

GiTiT T Grin T ; (14) 

Gt±(1',2';1,2) = -(raT±(l',2')4±(l,2) 



G 



TTTT.Ua 



(15) 



where we have dropped the argument (1', 2'; 1, 2) for ease 
of notation. Using ljll|l . we find that the spin and orbital 
symmetry of the pair of particles is directly refiected in 
the vertex F: 

Gs/t(1',2';1,2) = -G(1)G(2) [<5(1' - 1) ± <5(1' - 2)] (27r)3 
+ (*/n)G(l')G(2')r±(l', 2'; 1, 2)G(1)G(2), (16) 

where T denotes one of the triplet states To^±. Prom Eq. 
Ill(i|l we see that the vertex (and therefore the scattering 
amplitude) has either a bosonic (symmetric) or fermionic 
(antisymmetric) behavior depending on the spin state. 
Therefore, we calculate the unsymmetrized vertex T giv- 
ing the i-matrix t and the scattering amplitude /, and 
(anti) symmetrize the latter according to the spin state. 

Prom now on we consider the Green function in mo- 
mentum and frequency space. Taking into account con- 
servation of momentum and frequency, the vertex satis- 
fies the Bethe-Salpeter equatior-^^ 

r(p',p;P) 

i 



ft(27r) 



HP\p\P) 

dkA{k,p; P)G{ki)G{k2)T{p' , k- P), (17) 



illustrated in Fig. Efa). We have introduced the ir- 
reducible vertex A, the intermediate momenta fci 2 = 
P/2zL k given in terms of the relative momentum k, the 
frequency ui and the notation p = (p, lu). All the possible 
exchange diagrams should be included in A while avoid- 
ing double-counting physically equivalent diagrams. We 
first consider zero temperature ksT = 0, and discuss fi- 
nite T effects later. 

In a first stage, we use the random phase approxi- 
mation (RPA) for the many-electron backgrounc^; this 
yields the screened interaction 



Viq) 



Vc{q) 



l-Vciq)x'{qy 



(18) 



given in terms of the bubble susceptibility with the 
momentum transfer q = (q, Wg) = p' — p . The RPA re- 
quires a high density, which is controlled in 2D via the 
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Figure 2: (a) The two-particle Green function and the Bethe- 
Salpeter equation for the vertex V. We only show the direct 
diagrams (i.e., without the exchange ones), (b) Lowest or- 
der diagrams (Ai,2,3) contributing to the irreducible vertex 
A. The wavy lines denote the screened Coulomb interaction 
V , given in RPA by resuming the bubble diagrams. 



parameter = meo/(?i^-y/7m) ^ 1, where n is the elec- 
tronic sheet density. First, one can consider the static 
limi^ 



Q III 

X = 



i-e(q>2Mi/i- ( ^ 



(19) 

because the dependence of with uoq is smooth, and can 
be therefore neglected in the integration of intermediate 
fermionic lines — which, as we shall see below, selects 
only intermediate states k at the Fermi surface: = 
^ hujq = - = for Cl = 0. Note that the 
divergence in the ojq — integration of V{q) ~ ^/{^q — i^p) 
near the plasmon frequency LOp disappears because it is 
an odd function of w. Next, the g-dependent part of 
X°(g,a;g — 0) in 2D vanishes at the Fermi surface, as g ~ 
2kF sin |6'/2| < 2^^. This justifies the standard Thomas- 
Fermi screening 



V{q)^ 



(20) 



with the screening momentum kg = 2me^/h^ = kFrs^/2. 

Within RPA, the renormalized one-particle Green 
function is given bj*^: 



G{k) 



Z 



Wfe - - zlml](fc) 



(21) 



for small uJk and £,k- In 2D, for GaAs and ksT^^k 0, 
one has the renormalization factosS^ifi .Z = l — rs(l/2 — 
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Xj 1^)1^/2 ~ 0.62, the renormalized massiS m* = rfi\V — 
ln(l/rs)rs/7r] ~ 0.96m entering in and the broad- 
ening (inverse Hfetime)i Im E ~ ^^In^. The latter 
vanishes for particles near the Fermi surface (^^ 0), 
which corresponds to well-defined quasi-particle states. 
For simplicitjiSi, we set Z = 1, m* = to, and therefore 
approximate the renormalized Green function by the free 
propagator 



G(fc) ~ Go(fc) = 



1 



(22) 



We now consider the irreducible vertex A in lowest 
orders in F,as mV/'KtT? ^ mel/kp ^ ^ 1. The lowest- 
order diagrams are shown in Fig. IS^b); they are the single 
interaction line Ai, the vertex renormalization A2, and 
the crossed diagram A3: 



Ai = V{q), 



(23) 



A2 = n9)(^ / ^'"i^^'^i'^) [^(ki - Pi) + - ki) 

(24) 



A, 



dkiB(ki,g)F(ki+pi)F(ki+p2), (25) 



(2vr)2 

with q = p[— pi and Q = p'^ — P\- The function 

n(k + q) - n(k) 



i?(k,g) 



(26) 



arises from the frequency integration of the bubble 
diagram, and involves the Fermi occupations factors 
n(k) = 0(— ^fc) at /cbT = 0. We can estimate A2 ~ 
l/(g)T/(fcj^)/dkii?(ki,g)/(27r)2 = y(g)F(fcF)x"(g) ^ 
—V{q)ks/kF (apart from negligibly small integration re- 
gions), and similarly for A3. Hence^^, A2,3/Ai ~ Ts <C 1 
and we can keep only the direct interaction line A ~ Ai = 
V , which corresponds to the ladder approximation. The 
criterion for the validity of the ladder approximation is 
usually expressed^ by Afc_F ^ 1. This low density regime 
is nevertheless consistent with RPA, as can be seen, e.g., 
from the Born approximation 



t~Vc^ Xkp 



1 


mVcikp) 


2 




2^ 


ft2 







< 1. 



In summary, we need to solve the following Bethe- 
Salpeter equation: 

T{p'^p;P) = Vip'^p) 

dkvCk - p)Goik)Goik2)T{p' - k P), (28) 



where ki 2 



B. Energy integration and logarithmic factor 

We now follow the derivation of the Cooper 
instability^, including a discussion of the less standard 
case of particles that are not in the Cooper channel, i.e. 
with pi 7^ — P2. Before solving Eq. lf28ll to all orders, we 
first consider its second-order iteration, r*^^^ (p' — p; P) = 
V{p' - p) + i/h{2Tr f J d~kV{k - p)GCki )GCk2)V{p' - k) . 
The Wfe-integration of the Green functions yields'* 



N{ki,k2) 



-ujk)Go{k2 



2 



+ 6 - ^1 - + 2iTjN{ki, k2) ' 



(29) 



with the function N{ki,k2) := 1 — n{ki) — n{k2)- The 
P-frequency HQ has been set to ^1 + £,2 to retrieve the 
i-matrix, and rj can be set to in the denominator. We 
now consider the Cooper channel p = Pi = — P2,k = 
ki = — k2, which gives q — |p' — p| = 2psin|0/2|, and 
n(fci) = n(k2) = e(-6) =^ N{k) = sgn(Cfe). This yields 
a single discontinuity in the numerator when = 0, 
which coincides with the zero of the denominator {£_k = 
when considering a vanishing excitation energy for the 
incident particles: ^ = ^1^2 — 0. 

As a consequence, the main contribution to the energy 
integration comes from virtual states at the Fermi sur- 
face, i.e. £,k — 0. We set k = kp in V, and integrate only 
on D{k) = D(k, k). The dominant term comes from both 
sides around k = kp and yields the factoi4fi 



1 r°° 

v.= — dkkD(k) 



2Trh-^ Ep' 



(30) 



This logarithmic divergence (as C ^ 0) plays a crucial 
role here. It allows us to neglect the fc— dependence of V, 
as it selects only intermediate states at the Fermi surface 
(their contribution is logarithmically dominant as ^ ^ 0) . 
This is the main effect of the fermionic background on the 
scattering, apart from the screening. 

We now consider finite temperatures with n{k) = (1 -t- 



)-\ For C 



KT) 



47rft2 



ksT ^ Ep, we find 

e {kBT? 



log 

knT 



E? 
J^p 



ksT + 4 



(31) 



(27) For ^ > fcsT" we recover l(3(Hl . while for ^ <C A:_bT we find 



2TTh^ 



log 



ksT 
Ep ' 



(32) 



This logarithmic factor II32II is the 2D equivalent of the 
one found in the discussion of the Cooper instability for 
both phonon-mediated^ or Kohn-Luttingei*i^ supercon- 
ductivity — with the 2D density of state V2d —m/ {2Tifi^) 
instead of the 3D one. 
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Different cut-offs arise if the Cooper channel condi- 
tion (P = 0) is not strictly respected. For instance, ex- 
periments might require a small but finite angle 2a — 
^(Pii ~P2) ^ 1 between the incident particles to prevent 
misalignment, in which case we have a total momentum 
P — 2pi sina ~ 2fci?Q! if pi = p2 and \pi — kp\, P ^ kp- 
Alternatively, particles might be injected in perfect op- 
posite direction (a = 0) , but with a different energy (this 
can arise e.g. in case of hot electrons, see Sec II V i.e. 
Pi ¥^ P2, leading to P = \pi—p2\- In these cases, we find 
at T = 



N{k, (j)) = e{k - k+) - e(fc_ - k), 



(33) 



where k± = y kp — (Psin0) ± Pcoscj), and (p is the 
integration angle for k [(/) = Z(k, p) -I- 7r/2 if a 0; 
(j) — Z(k, p) if pi 7^ P2]. In the Hmit \pi^2 — kp\ <C P we 
find 



(34) 



where we recall that P = 2kpa or \pi — P2\- 

Therefore, the divergence (due to the single discontinu- 
ity of the two Fermi surfaces for particles in the Cooper 
channel) is cut-oShy inax{£_/ Ep, ksT/ Ep, P/kp}. Away 
from the Cooper channel (P ~ kp), or for large tem- 
peratures {ksT,^ ~ Ep) the logarithmic factor disap- 
pears; in that case the fact that only virtual states having 
the Fermi energy {k ~ kp) contribute to higher orders 
does not apply. When pi^2 — kp, we must make sure 
that the ingoing and outgoing states are available, i.e. 
Pi,2,p'i 2 > ^F- For 7r/2-scattering and a ^ 0, we have 
Pi,2 — Pi{cosa ± sina) ^pi < kp{l + a) [see Fig. ^c)]; 
then a — P/2kp < £,/Ep. Similarly for px ^ P2, we need 
P2 = Pi — P > kp ^ P/2kp < ^i/Ep. As a consequence, 
the cut-off is determined by either ^ or fc^T. 

Before proceeding with the solution of the Bethe- 
Salpeter Eq., we comment on the difference with the stan- 
dard scattering theory in vacuum, where the scattering 
t— matrix is given by the Lippmann-Schwinger equationSl 
t™'^ + VCp'^Vp'' for an incoming energy E. As the 
single-particle Green function G'™'^(k) = 1/{E—Ek) is an 
odd-function around the divergence at Ek = E, no diver- 
gence develops, and no selection of intermediate states k 
occurs. 

In the case of a Fermi sea, the Lippmann-Schwinger 
Eq. is replaced by the Bethe-Salpeter Eq. (which can 
be written symbolically T — V + VDT), and G™'^ corre- 
sponds to the factor D, Eq. (l29ll . This is seen by setting 
Ep = 0, giving N{ki,k2) 1 and thus D Gf'ik) 
for the relative momentum k — (ki — k2)/2. The dif- 
ference between the Fermi sea and the vacuum cases Hes 
in the numerator N{ki,k2), which is a direct effect of 
Fermi statistics and assigns a negative sign to the ex- 
change processes where the transitions occur via two 
states ki,k2 < kp below the Fermi surface, as compared 
to the direct processes via intermediate states above the 
Fermi surface, ki^2 > kp (see Fig. In the Cooper 
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Figure 3: Scattering with a fermionic background, (a) Di- 
rect virtual transition: the initial states pi,2 first go to the 
available intermediate state fei,2 > kp (1), and then go to the 
final states (2). This process is represented by a term 
c^^,Ckc\cy. (b) Exchange process corresponding to c^^,c\ckCp: 
intermediate states with k-1^2 < kp first fill the final states 
p'1,2, creating a particle-hole excitation (1); the holes are sub- 
sequently filled by the initial states pi,2 (2). The effect of 
the many-particle fermionic background manifests itself in 
the function N{ki,k2) in Eq. 12911 . which adds a negative 
sign for the exchange process, and is therefore responsible, in 
the Cooper channel, for the logarithmic term 18011 that se- 
lects states at the Fermi surface, fci,2 — kp. In vacuum, we 
have only case (a), so N{ki,k2) = 1 and all states contribute 
with the same sign; this yields no divergence, and therefore 
no selection. 



channel, this factor is responsible for the selection of 
intermediate states at the Fermi energy (fci,2 — kp), 
via the log-divergence ll371|l : the latter arises because 
N{k) = sgn{^k) impHes that D{^k) l/ICfe| is an even 
function (after setting ^ = 0). We emphasize that the 
selection of virtual states at the Fermi energy (fc ~ kp) 
disappears away from the Cooper channel; see Eq. (IMIl . 



C. Solution as a Fourier series 

We can repeat the integration over the frequency and 
energy described above at every order. This results in a 
new Bethe-Salpeter equation 

m = V{e) + / dcl)V{cj,)t{0 - cf,), (35) 

with the screened 2D Coulomb potential at the Fermi 
surface 



i;(0) 



2TTe' 



2kp sin 10/21 + ks 



(36) 



To solve this integral equation, we expand w in a Fourier 
series: 

«(0) = y\ "«e™^ , t,„ = / d0F(0)e-"'^ (37) 

as well as t((/)). The solution of the Bethe-Salpeter equa- 
tion is then simply given by 



^ 1 - VVn 



(38) 
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This expression for the Coulomb scattering i-matrix of 
two electrons in the Cooper channel (p2 = — Pi), in the 
presence of a Fermi sea, is the main result of the paper. 
We note that the procedure followed here is not valid for 
very small r-g ^ £,/ Ep, ksT /Ep, which will be addressed 
later in AnnlBl 



D. Fourier coefficients of «((^) 

For the Fourier coefficients of we integrate lj^^7|l 

in the complex plane with z = e"^/^. We find 



cos(m7) 

kp cos 7 ^-^ 2n + m 

oddm>l 



E 



with 



sm7 



2fc, 



(39) 



(40) 



For numerical estimates, it is more convenient to write 
the result as 



2el 



kp cos 7 



£n = In 



(41) 



tan cos(2n7) - | sin(2n7), (42) 



Hence, the large n dependence is polynomial, w„ ~ n~^, 
which reflects the fact that the potential w(0) is non- 
analytic. This is in contrast with the 3D case, where 
the coefficients of the spherical harmonics decomposition 
areiS (27reg/fc|,)Q; (1 + 7-^^2-1/4), and their decay is ex- 
ponential in I {Qi is the Legendre function of the second 
kind). 



S. Small Ts approximation. 

We now expand Ij45|l in small Vg, and find 

.„c.-^|ln(2n7) + - 



nw-f 



^ln(2n7). 
kp 



(47) 

This expression is not valid for very large n, as we ex- 
panded to lowest order in wy. It remains finite in the 
limit ^ 0, because e§ ^ r^. 



IV. SCATTERING LENGTH AND EPR PAIRS 

We now apply our result to a realistic GaAs 2DEG, and 
study the dependence of the scattering amplitude on the 
scattering angle, Vg and temperature, before discussing 
the production and detection of spin-entangled electron 
pairs. 



2n-l 



An ^2 J2 



cos [7(2n — m)] 



odd m>l 

For instance, for n = we have 



m 



vo = 



In 



kp cos 7 



tan 



(i) 



(43) 



(44) 



1. Integral approximation. 

For small the sum l(39jl is smooth and can be ap- 
proximated by an integral, giving 

Vn o=i |sin(2n7) [J - Si(2n7 + 7) 

kp cos J I 12 



■cos(2n7)Ci(27i7 + 7)|, (45) 



where Si and Ci are the sine and cosine integrals. This 
expression can also be obtained (for large n) with a Hnear 
expansion of the sine in «((/)) ~ 27re§/ {ks + kp\(l)\) before 
calculating the Fourier coefficients. For realistic param- 
eters it is very accurate already for n > 1. It yields the 
asymptotics for n 1 



^0 



1 



2kpj'^ cos 7 



(46) 



A. The different scattering lengths 

We define the scattering length for singlets and triplets 

Xs/T{0) = \f{O)±f{e~7r)f, (48) 

following Ijl6ll . We recall the scattering amplitude f{0) 
defined in 



(49) 



with the t— matrix given by Ij38ll . Unpolarized sources 
contain 1/4 of singlets and 3/4 of triplets [see Eq. Q], 
which yields the scattering length 



(50) 



We also define the scattering length A*^^^ obtained from 
the Born approximation with the amplitude 



/(i)(0) 



MO), 



(51) 



ft2^/2^ 

as well as the corresponding bare scattering lengths Ac 
and a|^\ obtained by replacing i(0) with tciO) and vc{0), 
given by Eq. and 10. We point out that <r = rs/\/2. 
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B. Total scattering length 

We now take typical parameters for a 2D GaAs electron 
gasSl, = 13.1, — 0.86, and a sheet density n = 
IQi^m-^, and assume £, < ksT = IQ'^Ef {T = 20 
mK). First, we estimate the magnitude of the scattering 
and calculate the total scattering length integrated over 
tt: 

Atot = de\{e) = ^A's"* + ^A^°* = 3.39 nm. (52) 
with 

A*/ = 7.92 nm, A^°* = 1.88 nm. (53) 

This is consistent with the ladder approximation, which 
requires Atot^F = 0.54 < 1. We now use the Born ap- 
proximation lf5T|l and write A^ot — — Aox . We find 
for the direct part 

r de\f{e)\' 

Jo 

= Ap-^^^^ < 1 - 2sin7tan7Arcth I J^— I i(54) 
cos 7 1 \ y 1 + sm 7 / I 

with the Fermi wavelength Xp = 2iT/kp, and we recall 
that sin 7 = The exchange term is 

AW = Re / ddf{6)f{e - tt) 
Jo 

= Xp- — 1 log(sin7+^^ I — tan7Arcth (0037) 1(55) 
cos 2j I \ sm7 / J 

which yields a|q|j^j = 11.0 nm. We see here that the 
Born approximation significantly overestimates the ex- 
act result (I52II . For small , we can further approximate 
A^|] ~ Xprs/Vi and A^ ~ (AFr2/2) log(V2/r,), which 

gives an even greater length, A^ot — 17.2 nm. Note the 
divergence in the Ts ^ Hmit due to the forward scat- 
tering {q — 0) divergence of the unscreened Coulomb 
potential. 

C. Angular dependence 

We compare in Fig. ^ the angular dependence of the 
different scattering lengths for unpolarized electrons. We 
first see that the main effect of the Fermi sea is to reduce 
significantly the scattering by one order of magnitude 
compared to the vacuum case. The large renormaliza- 
tion is related to the relatively large value of = 0.86 
(and the large screening ks ~ kp), which strongly reduces 
the forward scattering divergence of the bare scattering 
Ac, A^^^(6' ^ 0). Furthermore, we notice that the Born 
approximation A*^^) is not vaHd in the Cooper channel, 
as higher order terms reduce the scattering amplitude. 




Figure 4: Scattering length \{9) for a GaAs 2DEG with sheet 
density n = 4 • 10^^ m^^ and fcsT = Bf/100 = 2 K, for unpo- 
larized incident electrons. We compare the exact result 15011 
for a 2D Fermi gas to its Born approximation A'^', and to the 
bare scattering result (with no screening) Ac and A^'; see 
Eq. H51|l. The main effects of the many-body background is 
to remove the 9 — divergence by screening the Coulomb in- 
teraction, and to significantly reduce the Born approximation 
via the logarithmic corrections in Eq. IHSII . The inset shows 
the ratio of the exact scattering length \(0) to the Born ap- 
proximation result A'^'. 



The fact that higher terms contribute significantly, de- 
spite the weakness of V{q), is due to their logarithmic 
enhancement by the factor v. 

The angular dependence of the exact scattering length 
A is similar, but not identical to the Born approximation 
result X^^\ as shown in the inset of Fig. ^ Importantly, 
X{9) is a smooth, monotonic (for < n/2) function, so 
that the interference mechanism survives for the produc- 
tion of EPR pairs at 6 = 7r/2. 



D. dependence and Born approximation 

We show in Fig. © a plot of the scattering length as a 
function of the density n or = meQ/h^^/Tm (top axis; 
we keep eo constant), for the angle 6 — 7r/2 (hence only 
the singlet channel contributes). There is a strong de- 
pendence ~ of the scattering, which could be stud- 
ied experimentally by varying n via a top gate. This 
dependence also roughly applies to A*^^) and Ac, while 
Ag^ ~ V^/kp ^ 

The Born approximation t{6) ~ v{9) is reached when 
the logarithmic factor disappears (v 0) and does not 
enhance higher order terms; see Eq. (EU. This occurs at 
high temperatures ksT Ep, for hot electrons ^ ~ Ep, 
or for electrons that are not in the Cooper channel {P ~ 
kp). The Born approximation is also reached in the very 
small Ts ^ 0.01 limit (not shown), when Wn ^ 1. For 
even smaller one can neglect kg inV, which yields the 
Born approximation of the bare Coulomb potential, t ~ 
Vc ■ We note that the effect of the Fermi sea is intrinsic 
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1.5 1 0.86 0.7 0.6 r 

\ \ 1 \ ^ 




Figure 5: Scattering length Aat6 = 7r/2asa function of 
the density n (see the corresponding Ts = meo/ft^\/™ on 
the top axis). The inset shows the comparison with the Born 
approximation and the bare scattering. 

in our calculation (by restricting the intermediate states 
to the Fermi surface), which therefore cannot recover, in 
the Ts ^ limit, the exact result for the bare Coulomb 
interaction Q. 

E. Dependence on T, ^ and a. 

Two effects appear when one varies the temperature 
T, the excitation energy ^ of the incident electrons or 
the impact angle a = Z(pi,— p2)/2 (a finite \pi — P2\ 
plays the same role). The first one is a change in the 
factor f appearing in the denominator 1 — vvn of the 
i— matrix ijSHl- For finite a, one should in principle in- 
tegrate 1^(0) ~ (TO/27r?i^) log (2asin |0|) over the inter- 
mediate angle (j) = Z(k, p) in the Bethe-Salpeter Eq. 
t{(j)) = v{(j>) + (l/27r) / d<t)v{(j))v{(j))t{4)). However, the 
dependence of i^(0) is smooth (logarithmic) compared to 
the behavior of the Coulomb potential at the Fermi sur- 
face, w(0) '-^ 1 / {4> + r . Therefore, we neglect this 
dependence and set e.g. (j) — ""/S in i^(0), which gives a 
constant v ~ {m/2'n:h^) log (a). Hence we take 

"^2^^^°s("^""{^'^'"})- 

The effect of this dependence on the scattering length 
A at = 7r/2 is shown in Fig. El In panel (a) we fix 
C/Sf = 10-^, ksT/EF = 10-2 and vary a (we recall 
that the Fermi temperature is Ep/ks = 162 K). First A 
is constant when a < ksT j Ep, and then increases slowly 
(the horizontal scale is logarithmic) as v decreases; for 
a— >l=>z^— >Owe recover the Born approximation A^^^. 
The function A is exactly the same if one interchanges ^ 
and fcsT (b), or permute all the three parameters ^, fc^T 
and a (c-f). 

Secondly, we take into account the requirement that 
the outgoing states p'j^ 2 ^-^e not occupied (and hence 





Figure 6: "Observational" scattering length \F At 9 — 1^/2 as 
a function of the impact angle a (a,b), the excitation energy ^ 
(c,d), and temperature T (e,f) (left axis). The Fermi occupa- 
tion factor _F = [1 — m(Ci)][1 ~ ''^{^2)] (right axis) enforces the 
requirement that the outgoing states ^J_2 = ? ± sin(2a) 
are available for the outgoing scattering states. The depen- 
dence of the "bare" A comes from u l|56|l and is the same 
in all graphs. The arrows indicate the position of the fixed 
values, e.g. S,/Ef = 10~^ fesT/^F = 10"^ (a), £,/Ef = 
10-'^, kBT/Ep = 10"^ (b), etc. 



available for the outgoing scattering states) , by introduc- 
ing the factor F = [1 - n(ei)] [1 - n{^'^)] = n(-C; )n(-e^), 
plotted in Fig. El (right vertical axis). At 9 = Tr/2 we 
have p[ 2 — Pi (cos a ± sin a) 2 = C i sin(2Q;) ; see 
Fig^c); hence for a large a the final state £,2 < will be 
already occupied and the scattering into this channel will 
be prohibited. The transition across the Fermi surface al- 
ways occurs at the largest quantity [e.g., at ksT ^ ^ in 
panel (a)], and is sharp when temperature is negligible 
(b,d). Note that here we consider that the initial states 
Pi,2 are always filled (i.e., with occupation 1), being ei- 
ther injected from the QPC or thermally excited. 

In (a) the maximum value of F (when aEp ksT) 
is F — 1/4, because the final energies |<^J 2I ^bT lie 
within the temperature window: n(0) = 1/2. The same 
occurs in (c) for ^ < fcsT, while F = 1 for ^ > ksT. 
For negligible temperature (b,d) F saturates to 1 when 
aEp <C ^. We finally note that (for the factor F) the 
panel (e) corresponds to the opposite of panel (a), (f) is 
opposite of (c), and (b) is opposite of (d). 

We can now consider the combined effect of A and F 
by defining the "observational" scattering length AF, giv- 
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ing the scattering length as could be measured in a real 
experiment. As a function of a, it only reproduces F 
by showing a smooth (a) and sharp (b) step. In (c), it 
first increases slowly because of the smooth transition 
from -F = l/4^1in the region of constant A (small 
before following the logarithmic increase of A (large £,). 
In (d) the transition is sharp and starts at XF = 0. In 
(e) the transition is smooth, and the measurable length 
XF follows A but is reduced by a factor of 4. In (f) 
there is an interesting non-monotonic behavior in the 
region above kBT/Ep > a; however it requires an ex- 
tremely small a, e.g. a ~ 0.2°, not reachable in a realis- 
tic experiment. We note that the rightmost parts of the 
graphs (above 10"^) are only indicative, because they do 
not correspond to regime assumed in the derivation of A 
{kBT/EF,^/EF,a^ 1). 

The scattering length vanishes logarithmically, A ~ 
l/log^vV) when all ksT,^ and a ^ 0. It is rem- 
iniscent of the vanishing of the inverse lifetime of a sin- 
gle quasi-particle excitation scattering with other elec- 
trons below the Fermi surface, when its excitation en- 
ergy vanishes^. However, the two cases are completely 
different: the lifetime diverges because of phase-space 
constraints due to the Fermi statistics; in our case the 
scattering of two particles above the Fermi sea vanishes 
because of the renormalization due to the Fermi sea. 




Figure 7: Scattering length for the bare Coulomb interaction 
[given by Eq . lO], compared to the classical value given by 

Ag = |/c(f)P + \fc(0 - 7r)|2. (a) <; = 0.6 corresponding to 

n = 




Figure 8: Ratio R(6 56) of the number of triplets/singlets 
collected in the detectors, placed at an angle 9 and with an 
F. Quantum oscillations aperture <56l = 5, 10° . 



In addition to the destructive and constructive inter- 
ference at — 7r/2, quantum oscillations can be seen 
in the bare scattering (of singlets, triplets or unpolar- 
ized sources) as a consequence of the angle-dependent 
phase <^ In I sin6'/2| appearing in Eq. 0. The number 
of oscillations is roughly given by = r^/v^, as illus- 
trated in Fig. [3 We see that the oscillations are ab- 
sent for n = 4 • 10^^ m~^, and only appear at lower 
density. At 61 = 7r/2 the quantum amplitude for un- 
polarized sources is half of the classical one given by 
Ag = \fc{d)\^ + \fciO - because the triplet con- 

tribution vanishes. 

For the many-particle result, such quantum oscillations 
could arise from the small imaginary part appearing with 

the logarithm in iy{T) = (\n ^ + f «^ , as it yields 

an angle-dependent phase when summing up the Fourier 
series. However, the phase is of order 1/ hi{kBT / Ep) <C 
1, and the oscillations are not visible. 



G. Production of EPR pairs. 

We now consider the setup of Fig. E^b), with detec- 
tors placed at an angle 9 ~ 71/2. The triplet channel is 
non-zero because of the small aperture angle 259 of the 
detectors. The scattering lengths for the singlet /triplet 



channels into the detectors read 

Xs/t{9,69)^2 f d9'\f{9')±f{n-9')\\ (57) 
Je-se 

which we use to define the ratio 

Nt 3Ar 

between the number Nt/s of singlet /triplets collected in 
the detectors. Here we have allowed for the case where 
the average scattering angle 9 deviates from 7r/2. We 
show a plot of R{e, 59) in Fig. for 59 = 5° and 
10°. We find very low values, i?(90°,5°) = 0.183%, 
or i?(85°,5°) ~ i?(90°,10°) ~ 0.7%, which shows that 
the collision entangler is efficient as singlets are predom- 
inantly collected in the detectors, even for 9 = 80°. 
We note that such deviation from 7r/2 would occur ex- 
perimentally, because the electrons injected through the 
QPCs have an angular spread. This spread could, how- 
ever, be reduced by the use of a lense-shaped top gate im- 
plementing a refractive medium for the electron motion^. 

Expanding in 59 we find 

Xs,T{9,59)^ 259{f±U^)^ 
+ ¥9^ [(/' T fL? - (/ ± /cx)(/" ± m + 0{59^) (59) 



10 



with 



(60) 



and ' denotes d/d0. For 
f" we get 



90° 



5°, and neglecting 



R{90°,69) = 



/(V2) 



56"^ = 0.178% 



(61) 



which is close to the exact value R = 0.183% found 
above. Thus the Taylor approximation is accurate, and 
it is clear that the triplet contribution can be made arbi- 
trarily small by reducing the aperture 59. We note that 
our calculation gives a ratio |/'//l of the order unity for 
a wide range of parameter, kBT/Ep = 10^^ — 10^^° and 
Ts = 0.1 — 1. Using the Born approximation we find a 
significantly lower value {59 = 5°) 



i?(90°,( 



1 



4(r, + 1)' 



t59^ = 0.05% 



(62) 



that would be more advantageous for EPR production 
(see the discussion on hot electrons in Sec. II V I|l . 



escape from the scattering region into the drain contacts 
(Fig. nj . Finally, we find the total scattering probability 
of two electrons 



Ptot = ^ - 3.4% 
w 



(65) 



(i.e., not necessarily into the detector). 

An additional interesting topic is the noised of the 
detected current. Since this is outside the scope of the 
present work, we only give here heuristic arguments. As 
the scattering probability is very small, one can assume 
that subsequent pairs do not interact with each other. 
This implies that the zero-frequency noise induced by 
the scattering should be mainly given by the partition 
noise 

Situ = 0)^IPs{l~Ps)^ IPs, 

which becomes Poissonian for P5 <C 1. The pulsed injec- 
tion of the electron via e.g. the periodic lowering of the 
QPC barriers reduces the stochastic nature of the tunnel- 
ing through the QPC if the lowering is sufficiently fast. 
On the other hand, this periodic change should lead to a 
more complex noise behavior for finite frequencies^. 



Hot electrons 



H. Current 

We now estimate the singlet current collected in the de- 
tectors for a given input current /. We neglect the angu- 
lar dispersion of the incident electrons (due to diffraction 
on the edge of the QPC), by assuming that the electrons 
occupy the lowest transverse mode in the QPC, and that 
the remaining spread could, in principle, be compensated 
by the use of lensing effect^. This gives longitudinal plane 
waves (with wave vectors pi,2) having a transverse width 
w roughly given by half the width of the QPC. 

We first note that the scattering length for the singlet 
channel is small. As = 0.24 nm. Taking w = 100 nm, we 
find the probability 

Ps = 7— = 0.06% (63) 

4 w 

for the singlets to be scattered into the detectors. First 
we assume that the electrons are injected simultaneously 
from the reservoirs (which can be achieved by opening 
and closing both QPCs at the same time); this yields a 
singlet current of 

Is = PsI ^ 0.6 pA. (64) 

We have considered a given current of / = 1 nA, which 
corresponds to a frequency in the GHz range for the open- 
ing and closing of the QPCs. Otherwise, the electrons are 
injected at random times, given by a Poisson process with 
rate Win = e//. Then the probability of finding two elec- 
trons inside the scattering region (i.e., in state |pi,P2)) 
is roughly P12 = (Wout/Wm)^, where Wout is the rate of 



It is interesting to consider the case of hot electrons 
with larger excitation energies ^ (e.g., a fraction of the 
Fermi energy Ep of the scattering region), obtained by 
applying a dc bias voltage AV across the input QPCs. 
It can be problematic to have incident electrons with 
such a wide range of energy, as this allows a mismatch 
of the incident energies (^1 ^ ^2) after averaging over 
both incident energy ranges, which introduces uncertain- 
ties in the scattering angle (see Fig . To avoid this 
situation, one can raise the QPC heights such as to al- 
low only a very small range of electrons to go above the 
QPC barriei»2-. For hot electrons with ^ ~ Bp, the ex- 
act result moves toward the Born approximation; see 
Eq. ll62|l and Fig. 31 Hence the scattering length in- 
creases [because the logarithmic factor decreases ^ ~ 
log{£,/ Ep)], while the triplet/singlet ratio becomes more 
favorable (i.e. smaller). On the other hand, the scat- 
tering length becomes smaller for higher momentum [as 
V ~ l/(fc + ks)]. Taking e.g. eAV = 3meV ~ Ep/5S., 
one finds values that are more favorable than for cold 
electrons: the singlet length is doubled A^"'^-' = 0.56 nm, 
while the triplet /singlet ratio is halved TZ — 0.10%. Note 
that hotter electrons have a smaller lifetime because of 
the increased phase-space that is available for scatter- 
ing with electrons below the Fermi surface^. Estimates 
of the electron-electron scattering length le-e have been 
obtained?, for a GaAs 2DEG using imaging techniques via 
an SPM, in good agreement with theoretical predictions^. 
In our case, one has le-e — 1-2 /im, which is similar to 
the scale L ~ 1 /im of our envisioned setup. Hence, one 
can expect some reduction of the signal due to relaxation 
into the Fermi sea, roughly given by ^ g-^/'e-e ^ 
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J. Detection of entanglement 

An important question is to demonstrate that the col- 
lected electrons are indeed spin-entangled EPR pairs. We 
propose here three ways to answer this question experi- 
mentally. The first one is to refocus the scattered elec- 
trons into a beam splitter and carry out noise measure- 
ments in one outgoing lead; in this situation enhanced 
noise (bunching) is a signature of the desired singlet 
state, while zero noise corresponds to entangled or un- 
entangled triplets"^. However, this method would prob- 
ably require some bridges to avoid the source reservoirs. 
The second one is to carry out tests of violation of Bell 
inequalitj«S2iiiii^, by measuring single-spin projections 
via a single electron transistor coupled to a spin filtering 
device. The latter can be a quantum dot in the Coulomb 
blockade regime^ or a QPO^ in a strong in-plane mag- 
netic field. The third method consists in adding a p-i-n 
junction^Sii^, allowing the recombination of the entan- 
gled electrons with unentangled holes into photons; one 
should then carry out the test of Bell inequality with 
the photons, by measuring their entangled polarization 
modes. 

In addition, we mention that the interference mecha- 
nism responsible for the vanishing of the triplets at 7r/2 
could be demonstrated by polarizing the incoming elec- 
trons spin, which can be achieved by applying a large in- 
plane magnetic field to turn the QPCs into spin filters^, 
or by replacing them by quantum dots^. The current 
recorded at the detectors, which is proportional to the 
fraction of incoming singlets ps = (1 — V^)/2, should 
then rapidly decrease as the polarization V of the spins 
increase. 



K. Creation of localized, non-mobile entangle- 
ment 

We discuss here a way to produce static spin-entangled 
electrons, described in Fig El We propose to replace the 
two "injection" QPCs by two quantum dots, each with 
an even number of electrons, so that two excess electrons 
are in the singlet ground stated. Lowering the tunneling 
barrier defining the dots allows for the simultaneous in- 
jection of one electron (of the singlet pair in each dot) into 
the scattering region. If these electrons are detected at 
7r/2 by the detector, we know that they must be in a spin- 
singlet state (with certainty 1 — Nt/Ns). As spin is con- 
served during the Coulomb scattering, the total spin for 
the two electrons left in the two dots must be zero, which 
corresponds to the spin singlet state. Therefore, one has 
created a locaHzed (non-mobile) entangled pair of elec- 
trons separated by an interdot distance I ~ 1 /xm. This 
instantaneous "creation" (obtained by post-selection) is a 
dramatic illustration of Einstein's "spooky action at a dis- 
tance". Contrary to the standard EPR paradox, it does 
not manifest itself in the results of measurement corre- 
lations, but in the "creation" of a non-local (in the sense 
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Figure 9: Creation of non-mobile entanglement. Each input 
QPC is replaced by a quantum dot {L and R) containing two 
excess electrons. The ground state in each dot is the singlet 
with total spin Sl -I- S^^/ = 8^-1- Sni = (a). One injects 
one electron from each dot (e.g. L' and R'), and allows them 
to scatter. When they are detected at a 7v/2 scattering angle, 
we know that they are in the singlet state Sj^/ -f Sjj/ = 0. (b). 
As the total spin is conserved, the two remaining electrons 
in dots L and R are also in the singlet state Sl + Sr — 0, 
and therefore form a localized EPR pair, whose members are 
separated by the interdot distance 

of non-overlapping wavefunctions) quantum state. Such 
process, similar in some sense to entanglement swap- 
ping or quantum teleportation, could be useful in a scal- 
able quantum computer to create entangled pairs without 
having to go through the standard sequence of swapping 
state, which requires to move one electron to neighboring 
site of the other one, entangle them by local interactioni^ 
and move the electron back to its original place. How- 
ever, we note that our scheme requires a large number of 
collisions, of the order of l/Pg; this number scales fortu- 
nately more slowly (~ 1/^^) with the aperture 60 than 
the precision (- Ns/Nt ^ l/SO"^). 



V. ELECTRON-PHONON AND ELECTRON- 
ELECTRON INTERACTION 

In this section we investigate the question whether 
phonons can influence the scattering amplitude in a sig- 
nificant way. We first note that the scattering of electrons 
on real phonons can be neglected here, as it is strongly 
suppressed at low temperature. This is illustrated for 
instance by the absence of phonon effects in the experi- 
ments of Refiii^. However, the effective electron-electron 
interaction arising from the exchange of virtual phonons 
does not depend on temperature, so that it could play a 
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role in the electron-electron scattering. Our goal here is 
to estimate it and compare it to the screened Coulomb 
interaction which we have considered so far. We shall 
see that the contribution of acoustic phonons (deforma- 
tion and piezoelectric coupling) is negligible, while the 
polar phonons give a smooth monotonic decrease of the 
electron-electron interaction that is less than 20%, and as 
such does not change qualitatively the results presented 
in Sec. Id 



A. 2D phonon-mediated electron-electron inter- 
action 

In 3D, the effective electron-electron interaction is 
given bji^ 



ph 



2V ^ ^ k+q k'-q 
k,k' q 



A. 



(66) 



where V the normalization volume, q — {q^,qz),k,k' are 
3D vectors and q, k, k' are 2D vectors in the plane of the 
2DEG (in the following we keep the notation q — |q|). 
The electron-electron interaction matrix element reads 



WP'^iq) = |M(g)p 



t^ph(g) 



hV LU^-iO^iq)' 



(67) 



where ujph{q) is the phonon dispersion, and M{q) is the 
matrix element for the electron-phonon interaction 



(68) 



k(7 



Here and ct are phonon and electron creation opera- 

q k 

tors. We shall consider the lowest order in W^^, which 
in the Cooper channel allows us to take the static limit 
w = as all the energies involved in the scattering are 
the same, Ei ~ Ep- 

The electron-phonon interaction M{q) as well as the 
effective interaction WP^{q) are always 3D as they in- 
volves coupling of the 2D electrons with the bulk 3D 
phonons. There is no 3D screening of the bare ion- 
electron Coulomb interaction, as there are no mobile 
charges in the bulk. Now we define an effective 2D inter- 
action W2Li(q), which we shall compare to the unscreened 
2D Coulomb interaction Vc- We assume that the elec- 
tron wave function is separable into a plane wave |k) and 
a confined lateral function For instance, one can take 
an infinite square well of width L 



i'iz) = I; sin 
which yields the width of the 2DEG 



1/2 



L 



1 

2^ 



(69) 



0.18L. (70) 



We prefer to consider the alternative variational solution 
of the triangular well present at the interface^. 



with the width 



|^(z)p ^ '-n^z^e-^^ 



d = , 

K 



(71) 



(72) 



as it allows for simple analytical expressions. We define 
the effective W^^iq = k - k') by 



and gefei^ 



'2D\ 

(73) 



W^lS(q) = JdzJ dz'\^P{z)\^mz')\^W^\q;z-z') 

= ^ j dq.W^\ci,qz)\I{qz)\\ (74) 



with WP*^(q,z) = (l/27r)/dgW^Pi^(q,g^)e-'9^ and the 
form factor 



I{q,) = / dze'9^^|V(z)|^ 



(75) 



The latter is particularly simple for the triangular well, 
I{lz) — {iqz/i^ — 1)"'^- Our goal is to find the strength of 
this additional e-e interaction relative to the unscreened 
Coulomb potential Vc, by defining the ratio 



ph 



2D 



Vc 



[q = kf 



(76) 



Parameters for GaAs We consider a well of width 
d — 5 nm, and take the following parameters^ for GaAs: 
the mass density p„i — 5320 kg/m'^, the deformation 
potential constant D = — 7eV, the piezoelectric con- 
stant ehi4 = 1.44 • lO^eV/m, the acoustic sound veloc- 
ity c = 3700 m/s (we assume here that c is the same 
for both longitudinal and transverse phonons), the op- 
tical (longitudinal and transverse) phonon frequencies 
hujLO — 36.6 meV, huixo — 33.8 meV, the ionic plas- 
mon frequency: flpA — 85.5 meV and, finally, the low- 
and high-frequency dielectric constants e(0) =12.9, e(oo) 
=10.89. 



B. Acoustic phonons: coupling to the deformation 
potential 

We first consider electrons coupled to the acoustic 
phonons via the deformation potential. The electron- 
phonon matrix element is^ 



(77) 
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where D is the deformation constant, pi is the mass den- 
sity, and the dispersion relation is = c\q\. The 
static effective e-e interaction is a constant 



which yields in 2D for the triangular well iTtTIi 



<D(q) 



D2 3k 

PiC^ 16 



The ratio Ij7(i|l becomes 

r{q) = -q 



:^i^«=i-_i.4.io-3 



(78) 



(79) 



(80) 



which shows that the effective interaction T^j^d can be 
neglected for deformation potential coupling. 
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Figure 10: Effective 2D e-e interaction Wj'd from LO polar 
phonons as a function of the scattering angle 9. We take a 
well with a width d — 5nm, and fix g = 2fcFsin(S/2) and 
LO — Q. We compare it to the unscreened Coulomb interaction 
Vc- Inset: Ratio TZ = WJd/^c from 19111 ; the horizontal line 
corresponds to the small-g approximation IHOII . 



C. Acoustic phonon: piezoelectric coupling 

For piezoelectric coupling, the matrix element readsSS 



(81) 



with the polarization constant e/ii4 and the anisotropy 
factor 



m 



9q2g4/2|g-l6 (LA) 
+ g^)/4kT (TA) 



(82) 



for longitudinal (LA) or transverse (TA) phonons. It can 
be replaced by Ala = 0, ^ta = 1/4 for a 2D system 
constraining momentum transfers to qz — 0. This gives 



%Pic\q\ ' 



(83) 



The static e-e interaction is therefore proportional to the 
3D Coulomb interaction: 



W'P'^iq) 



1 

Pi \2cerJ q^ 



(84) 



Performing the transformation Ij74|l . we find for small 
q < kp/\Q the effective 2D potential: 



w^jS(q) 



1 / hii \ e 



Pi \2cer J 2q 



2 ^2 



(85) 



This corresponds to the 3D 2D transformation of the 
Coulomb potential, i.e. 1/|<7|^ 2/q . Finally we get 

r = -^(^y = -1.5-10-^ (86) 
Pier \ c J 

Hence the 2D piezoelectric contribution is also negligible. 



D. Optical phonons: polar coupling 



The electron-phonon matrix element is^ 



M{q) 



e(oo) 



1 



27rp2 

-j^hcoLoV, (87) 



where e(0) and e(oo) are the static and high-frequency 
dielectric constants. For uj — 0, this yields 



WP^iq) ^ - 



1 



e(oo) 



1 



Hence the effective 2D potential is, for small q, 



W!^{q)^^ 



and we get the ratio 



1 



e{oo) 



1 



27reg 



e(oo) 



1 



= -19%. 



(89) 



(90) 



As the ratio \r\ is rather large, it is important to con- 
sider here the more accurate expression, valid for larger 
q, found by performing the integration (I74|l with the tri- 
angular well solution (jjlj: 



e{oo) 



1 



27reg 



\{K^~q^y 



{-3q^ + lOq^K^ - Ibqn^ 



(91) 



which is plotted in Fig. EH We find the rather unex- 
pected result that virtual (optical) phonons give a signif- 
icant contribution: the ratio is \r\ < 20%. However, the 
effect of is monotonic and will not change qualita- 
tively the scattering of two electrons in a GaAs 2DEG. 
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VI. KOHN-LUTTINGER INSTABILITY 

Having found the scattering vertex in lowest order, 
we now consider higher-order diagrams and examine 
whether superconducting fluctuations could have an ef- 
fect on the scattering. It has been known for a long 
timei^ that in 3D the second-order crossed diagram A3 
in the irreducible vertex (see FiglJl can lead to a pairing 
instability and a transition to superconductivity. The 
origin lies in the susceptibility x° entering A2,3; being 
non-analytic, its spherical harmonics have a polynomial 
asymptotic decay x? ^ with respect to the coefficient 
I of the spherical harmonics decomposition, while the sin- 
gle interaction is analytic and therefore yields vi ~ e~K 
As x° oscillates, the irreducible vertex becomes attractive 



for sufficiently large I: A; 



-Aa^i < 0. The transition 



temperature is found from the Cooper divergence of F, 
i.e. by the relation 1 = v^dAi, where u^u is the same 
as ll32|l with the 2D density of state m/2'Kfi? replaced by 
the 3D one. This yield an infinitesimal temperaturei^, 
ksTc ~ exp(— 10^) for a metal with = 4.5. 

In 2D the equivalent transition does not occur, be- 
cause there is no instability for particles below the Fermi 
surface: q <2kp and x° = —-^0 has no (negative) har- 
monics. It has been shown, however, that higher-order 
diagrams (in A) can lead to a transitioni^. Alternatively, 
finite energy transfers can induce "pseudo-pairing" in the 
rf-wave Cooper channeU^. In our work, however, we are 
not interested in a superconducting transition; rather, we 
would Hke to verify that the scattering vertex for the in- 
jected particles (which are above the Fermi surface) is 
not substantially renormalized by the standard (lowest- 
order) Kohn-Luttinger instability with no energy trans- 
fer. 

The singular part of A2,3 originates from the function 
B{k.,q) evaluated near g = 0. Neglecting the variations 



of V in the g-integrals (I24I25|I (this corresponds to ap- 
proximating V{r) by a very short range potential, e.g. a 
(^-function) , we evaluate V at the singular points of B 
and write 

A2 ~ 2V{Q)V{2kF) j dkiB(ki,g) - -V{0)V{2kF)x\q) 

(92) 

A3 ~ V\Q) j dkii?(ki, Q) = ~V^{Q)x''{Q)/2 (93) 

where q = |p' — p| ~ 2kF and Q = |p' -I- p| — 2fc_F. We 
see that we can neglect A2 ^ r^Aa; we also take the static 
limit itToll of x^{^) cis it varies on a scale ~ Ep- As Q 
should be sHghtly above 2kF, we take 



Q = kF(2 + c) cose/2 



(94) 



with c = (,/Ef 2{p - kF)/kF < 1 and 16*1 < 1. The 
Fourier coefficients of the crossed diagrams are therefore 
given by 



(l + c/2)2cos2(6'/2)' 



(95) 



with 9 ~ v2c. In lowest order in c, we find for the 
singular part 



A3,„--^'(0) 



27r;i2 4^ 



and the asymptotics 



A 



3,n 




— sin(n6'- 7r/4), 



(96) 



(97) 



comparedi^ to l'"^ in 3D. Note the oscillatory behavior, 
which allows for negative values. The instability temper- 
ature can be estimated by requiring 



|A3,n| >Vn^n> no 



(98) 



where we neglected the oscillating sine function, and used 
lUfijl for Vn. We find 



/cbT ^ Efc 



-1/|A, 



3. .10 



(99) 



Note that the parameter c appears in 2D because of the 6 
function in (jl9ll . and is absent in 3D. For GaAs and tak- 
ing c = 0.02, we find ksT/EF ~ 6"^"°, which means that 
the attractive effect of the crossed diagram is completely 
neghgible, and does not lead to any sizable fiuctuations of 
the scattering vertex. For a metal with Ts ~ 4.5, the tran- 
sition temperature in 3D was foundli to be ^ Epe^'^^'^'^^; 
in 2D we cannot neglect the sine as in l(98jl : numeri- 
cally we find the temperature ~ Epe'^^^'^. This can 
be larger than in 3D for c > lO"**, despite the fact that 
the asymptotic decay of Vn ~ n~2 is much slower than in 
3D {vi - e-'). 



VII. CONCLUSION 

The prospect of experiments probing individual elec- 
tron collisions in a 2DEG is a natural motivation to study 
the problem of two electron interacting via Coulomb in- 
teraction in the presence of a Fermi sea. One of the main 
result of this work is the expression II38|I for the scatter- 
ing amplitude for two electrons in the Cooper channel. 
We found that the presence of the Fermi sea yields a sig- 
nificant renormalization of the strength of the scattering, 
rather similar to the renormalization found in the discus- 
sion of the Cooper instability. This is closely linked to 
the selection of intermediate states at the Fermi surface. 
Away from the Cooper channel, this selection disappears, 
and the Born approximation is valid. The overall angu- 
lar dependence is fairly unmodified and smooth. There 
is a sizable dependence on the sheet density, while the 
dependence on temperature, energy and impact angle is 
strongly influenced by the Fermi occupation factors. The 
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total scattering length is Atot — 3nm, which is of the 
same order as the Fermi wavelength. 

We discussed how to use such collisions to produce 
EPR pairs at a scattering angle of 6* = 7r/2. This mech- 
anism is rather robust against imprecisions in 0, for an 
output singlet current around 0.5 pA. The EPR produc- 
tion was found to be slightly more efficient in the case 
of hot electrons. We discussed detection of entanglement 
and quantum interference, and proposed a way to create 
localized EPR pairs separated by mesoscopic distances. 

We studied phonon-mediated electron-electron inter- 
action. We found that the dominant contribution comes 
from polar coupling to optical phonons, but does not af- 
fect qualitatively the Coulomb scattering. The strength 
of the Kohn-Luttinger superconducting instability was 
calculated, and shown to be negligible. Finally, we devel- 
oped (in AppendixEl an alternative calculation valid for 
diverging forward scattering contributions, and showed 
them to be negligible in GaAs. 
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Appendix A: SMALL r,-APPROXIMATION AT 

= 7r/2 



and neglect the second term which is order of . Then 



i(7r/2) 



2v 



(1 + lyvlogrg)'- 



?i>0 



1 



1 



and finally 



tin/2) 



2v 



(1 + lyvlogTs)"^ 2V2' 



4n + l 4n + 3 ^ 
(A5) 

(A6) 



This approximation is good for very low r^; the error is 
~ 10% for rg < 0.09, which corresponds however to a 
very high density n = 4 • 10^^. 

We proceed similarly for the derivative i'(7r/2): 



t'{e/2) = 2 {-ti + 3t3 - 5t5 + 7t7 -9tg + . 



(A7) 



4w 



1 



1 



An-1 4n -I- 1 



(A8) 

The sum yields -7t/8V2+ (1/4)(-1)^ with N 00. 
Neglecting the oscillating termS?, we get 



t'in/2) - -t(7r/2). 



(A9) 



We find in this approximation a very simple form for the 
ratio 



Ri0, 



1 



(AlO) 



Here we derive analytical expressions for the i-matrix 
and its derivative at 6* = •7r/2, which we then use to com- 
pute the ratio R. We first note that t{-n/2) is an alter- 
nating series 



t(7r/2) =to + 2^(-l)"i2n 



(Al) 



n>l 



where i„ = Vn/{'^ — vvn)- We write the differences as 

tn-tn+2 = hn/[{l-VVn){l-VVn+2)\ with /l„ = Vn-Vn+2, 

which allows us to get the smoother series 



t(7r/2) 



{I - VV(^){1 - VV2) {l-VV2){l-VVi) 



(A2) 

Now we use for the denominators the very small Ts ap- 
proximation (j47ll Vn ~ — wlogr^, V = 2e§/fci?, giving 



t{lT/2) = 



1 



(1 + VV logTsY 

For the numerators we write 

, cos 7 cos 37 

hn = 2v 
+2w 



{ho-h2 + hi-hQ + ...). (A3) 



2n+l 2n + 3^ 
cos 57 — cos 7 cos 77 — cos 37 



2n + 5 



2n + 7 



(A4) 



This corresponds to the Born approximation result, Eq. 
II62II , in the limit of no screening (ts 0) . This is some- 
what surprising, as our result ijAfill still contains both the 
screening (finite Vg) and the resumed higher order terms 
(responsible for the term vv in the numerator). One must 
further expand i(7r/2) in small in order to recover the 
Born approximation with an unscreened potential 

t{Tr/2)c^v^^rsn—^Vc{q^kFV2). (All) 
V 2 

The Fourier series is not well defined in this case, because 
of the forward-scattering divergence Vc(0). We also note 
that the number of Fourier coefficient required to reach 
convergence of the numerical Fourier sum increases dra- 
matically to rimax — 70000 for n = 10^'^ (with the heuris- 



tic dependence n„ 
more peaked. 



-0.8^ 



as the potential becomes 



Appendix B: FORWARD SCATTERING WHEN 

Here we consider carefully the limit of vanishing 
0, by following a different approach to solve the Bethe- 
Salpeter Eq. that allow us to study the contribution of 
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forward scattering states. These are indeed important in 
the very small limit, as the unscreened Coulomb scat- 
tering cross-section has a forward scattering divergence 
(i.e., for vanishing momentum transfers, g = 0) in 3D 
and 2D. 

The calculation which was presented in Sec. IIIII is 
based on the logarithmically dominant contribution of 
u ~ logc with c = max{kBT,^)/Ep; it is a many- 
body effect, related to the sharp edge of the Fermi sur- 
face, that occurs only in the Cooper channel P2 — — Pi- 
This approach fails for in the situation where is very 
small, when the screening is too small to reduce the for- 
ward divergence of the unscreened Coulomb potential. 
In this situation, one must consider carefully the con- 
tribution of forward scattering intermediate states with 
g=|k — p|~Oas they yield at large term V ~ 2TTe'^/ks. 
For such states, one must keep the restriction k ~p ~ kp, 
but consider the contribution of small angles = Z(k, p), 
defined bj*^ 



d 
p 



5 
P 



P 



kf 



« 1, (Bl) 



where d = k — kp, 5 = p — kp and 9 = Z(p,p') is 
the scattering angle. This corresponds to the forward 
scattering into virtual states; we do not yet specify the 
real scattering angle 9 between the initial and final states. 
We now go back to the iterations of the Bethe-Salpeter, 
after the frequency integration , still considering the 
Cooper channel only. 



1. Angular integral. 



We introduce the polar notation = {ki^ipi) and write 
the potential V{ki,(l)i;k2,(j>2) = VCki — k2), assuming 
p = (p, 0),p' = {p,9). We first consider the angular 
integrals for the n"^ order iteration of the Bethe-Salpeter 
equation: 



1 

2^ 



d(f)iV{p,9;ki,(j)i)— / (i(/)2V^(fci,(^i;/c2>2) X...X 



2tt 



2tt 



d4>n-lV {kn-2, 4'n-2\ kn-1, (j^n-lW {kn-l, (j^n-l] P, 0) 

(B2) 



V{p, 9- fci, 0)y(fci, 0; fc2, 0)... F(fc„_2, 0; fc„_i, 0)T/(fc„_i, 0; p, 0) 



n— 3 



= fe(fc„_2, fcn-l) 



2. Energy integral. 



First we consider 



For the integration of a{k) in A := 
{I / 2Tr) J^^'^^^ dkkD{k)a{k), we set k' — p and ex- 
pand the result; this gives 



b{k',k) 



2^el 



27reg 



k\ fc.. 



-P\ 



(B3) 



We truncate the energy integration to a small range 
k € [fci? = p — (5, p + (5] around p and define B{k') := 

(1/27:) JI:^^^^ dkkD{k)b{k',k). We can expand the 
lengthy result in lowest order in ks and find 



9 4 



h?kJk' 



■P\ 



-[l + 0{kjp)]. 



(B4) 



Thus, the subsequent fc'-integration will be mainly given 
by a small region around k' = p, as expected. Hence 
one should set fc' = p in the exact result B{k') before 
expanding in lowest order of ks , which yield 



B{k' ^p)~i- 



h'^k^ 



[1 + OikJp)]. 



(B5) 



A 



^^^[l + 0{ks/p)]. 



2h^k, 



(B6) 



For the last term containing the scattering angle, Ae := 



J^Ik^ dkkD{k)ag{k) we have 
A 



,Ji^ + OikJp)]^A^ (B7) 



with q = 2fcFsin|6'/2|. Writing the series t = E„ 
we find for the rt*''-order term (n > 0) 



27re, 



nn-l 



'-D 



with 



D = — c= ^. 

Stt Stt Ep 



(B8) 



(B9) 
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3. Result 

So far we have only considered the case where all 
intermediate angles are 0^ = 0, i = 1, ..n, and only 
the final angle is the scattering angle 9. There are 
n equivalent cases giving the same contribution, where 
the first j angles are zero and the remaining ones are 
(pi = 9,j < i < n. Thus the expansion series of the 
Bethe-Salpeter Eq. is t = ^ + Z]n>o Perform- 
ing the summation, we get the scattering i— matrix for 
Coulomb scattering in the limit of very small <C 1: 



t = 



2-D fe,«g 2TTel 1 



[l-DY 



q [l-Dy 



(BIO) 

with q = 2kp sin \ 0/2\. We can further expand this result 
in -D <C 1 t — 27reg/q — Vc; however, this merely 
refiects the fact that the Born approximation with un- 
screened potential is accurate in the limit ^ be- 
cause eg is also proportional to r-g. We also note that 
with this result the scattering amplitude has a phase 
2Arctan|Z?| = 2Arctan((^/87r); however, it is independent 
from the angle 9 and therefore does not yield quantum 
oscillations in the singlet /triplet scattering length as dis- 
cussed in Sec. IIV Fl 

In the case of forward scatteringSi with 9 = 0, the 
counting of equivalent arrangements of the intermediate 
states gives a factor 2" instead of n. The result has a 
more familiar form: 



t ■■ 



27reg 



1 



1 - 2D 



(Bll) 



As expected, the forward scattering amplitude ijBllll is 
larger than ijBlflll by a factor ^ q/ks ^ l/rg. 



4. Comparison 



We now compare the n'^ order obtained here, ti^^ de- 
fined in Eq. ljB8|l with the calculation of Sec. IIIIL i.e. the 

m}^ Fourier coefficient tm'' = Vm{vvm)"'~'^ ■ Their ratio 
is 



in) 



^(") 



n27reo 



1 



log 



kp 



kg logc 



kp c 

(B12) 

where have considered the small Vg approximation ijTzjl 
Vm — (27reo/fcF) log (/cs/^f). Therefore, the calculation 
of Sec. mil is valid provided that 



Ts I log Ts I > 



V2 |logc| 



(B13) 



Taking c = S/kp = 0.01, we find that we need > 
0.01, which is always the case for typical semiconductor 
material. In conclusion, the large value of Ts does not 
allow one to see the contribution of forward scattering 
into (virtual or real) states. 
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